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ABSTRACT 

Aims. We describe here a new full 2-D parameterization of the PSFs of the three XMM-Newton EPIC telescopes as a function of instrument, 
CIh energy, off-axis angle and azimuthal angle, covering the whole field-of-view of the three EPIC detectors. It models the general PSF envelopes, 
Q the primary and secondary spokes, their radial dependencies, and the large-scale azimuthal variations. 

» I , Methods. This PSF model has been constructed via the stacking and centering of a large number of bright, but not significantly piled-up point 
' sources from the full field-of-view of each EPIC detector, and azimuthally filtering the resultant PSF envelopes to form the spoke structures 
' and the gross azimuthal shapes observed. 
' ' Results. This PSF model is available for use within the XMM-Newton Science Analysis System via the usage of Current Calibration Files 
XRT/_XPSF_001 l.CCF and later versions. Initial source-searching tests showed substantial reductions in the numbers of spurious sources being 
detected in the wings of bright point sources. Furthermore, we have uncovered a systematic error in the previous PSF system, affecting the 
entire mission to date, whereby returned source RA and Dec values are seen to vary sinusoidally about the true position (amplitudewO.8") with 
' source azimuthal position. 

Conclusions. The new PSF system is now available and is seen as a major improvement with regard to the detection of spurious sources. The 
new PSF also largely removes the discovered astrometry error and is seen to improve the positional accuracy of EPIC. The modular nature of 
0^ the PSF system allows for further refinements in the future. 

o 
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• ^ 1 . Introduction this shared equally between the paraboloid and the hyperboloid 

r\ 1,. /. , ■ ■ halves of the Wolter I configuration. The maximum diameter 

XMM-Newton (Jansen et al. 2001), a cornerstone mission of . , , , • • ■ „ ■ • ■ 

13 . \,„„„ . , • , mirror shell IS 70 cm, and the outer and inner mirror shell thick- 

ESAs Horizon 2000 science program, was designed as an 

X-ray observatory able to study cosmic X-ray sources spec- 
troscopically with the highest possible collecting area in the 
0.2-lOkeV band. This high throughput is achieved primarily 
through the use of 3 highly nested Wolter type I imaging tele- 
scopes. The design of the optics was driven by the need to have 
the highest possible effective area up to lOkeV, and in partic 
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nesses are 1 .07 mm and 0.47 mm respectively (Aschenbach et 
al. 2000). 

One of the three co-aligned XMM-Newton X-ray tele- 
scopes has an unimpeded light path to the primary focus, 
where the European Photon Imaging Camera (EPIC) pn cam- 
era (Striider et al. 2001) is positioned. The two other tele- 



1 -71 1 1 1- . • 11 • -jz scopes have Reflection Grating Assemblies (RGAs) in their 

ular at ~7 keV, where the K lines of astrophysically significant , , ,.rr ■ r , ■ ■ ,■ • 

light paths, diffracting part of the incoming radiation onto their 



iron appear In grazing incidence optics, the effective area is 
generally increased by nesting many mirrors together and pack- 
ing the front aperture as much as possible. In the case of the 
XMM-Newton miiTors, each of the three telescopes contains 
a mirror module of focal length 7.5 m, comprising 58 nested 
mirror shells, the axial length of the total miiTor being 60 cm. 



secondary foci (where the Reflection Grating Spectrometers 
(RGS; den Herder et al. 2001) are situated), leaving the remain- 
der to travel straight through to the primary foci, where the two 
EPIC-MOS cameras (Turner et al. 2001) are positioned. 

A critical parameter determining the quality of an X- 
ray mirror module is its ability to focus photons, i.e. its 
Send offprint requests to: A.M. "Rssid Point Spread Function (PSF). Each of the three Wolter type 
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I X-ray telescopes on board XMM-Newton has its own 
PSF. As an example, FiglT] shows the in orbit on-axis PSFs 
of the MOSl, MOS2 and pn X-ray telescopes, registered 
on the same on-axis non-piled-up source (specifically of 
2XMM J130022.1H-282402 from ObsID 0204040101 from or- 
bital revolution 823). This figure shows the shape of the PSFs, 
with the characteristic radial spokes. Also note the coarser 
larger-scale shapes, in particular the strong triangular form of 
the MOS2 PSF, and the weaker pentagonal form of the MOSl 
PSF. Note also that these coarse shapes are quite diff'erent in 
the 3 EPIC cameras. 

Any measurements of the PSF by the EPIC cameras may 
depend on the instrument readout mode, through combinations 
of out-of-time (OOT) event smearing and/or pile-up. The PSF 
can be severely affected by pile-up effects when the count rate 
exceeds a few counts per frame. Depending on the selection 
of event types in the EPIC event analysis process, a hole can 
even appear in the core of the PSF due to photon events under- 
going severe pattern migration into unrecognized and rejected 
patterns, or yielding a reconstructed energy above the on-board 
high-energy rejection threshold. This is seen in Fig|2] (right), 
where the PSF of a severely piled-up M0S2 on-axis source is 
shown (specifically of GX 339-4, ObsID 0204730301, revolu- 
tion 783). 

Much more outer detail can be seen in the piled-up PSF 
(see FiglJl right). Note for instance how the piled-up central 
hole follows the general MOS2 triangular structure. The 16 
primary radial spokes are caused by the spider structure (see 
Fig|2l left) supporting the mirror shells (the primary spokes 
in the image actually lie between the spider 'legs', i.e. were 
the spider absent, the PSF would be as bright as the primary 
spokes around the whole azimuth). Note that the spider support 
structure at the paraboloid front aperture is the only support 
structure for the entire module; there is no equivalent coun- 
terpart behind the hyperboloid rear aperture. Looking further, 
note the 16 lower-intensity secondary spokes lying between the 
primary spokes. These are thought to be due to low-level scat- 
tering from the sides of the spider legs. The dark lanes visible 
at larger radii from the secondary spokes are due to the electron 
deflector which is mounted after the rear aperture of the mirrors 
and whose legs are aligned with those of the front-end spider. 
The coarser (triangular and pentagonal) image structures are 
believed due to deformations in the mirror shells, certain sets 
of mirror shells, believed to be the outer shells (de Chambure 
et al. 1999), not being perfectly circular in the 3 EPIC mirror 
modules. Finally, though not formally due to the mirror sys- 
tem, note the diffuse streaks of OOTs to the top and bottom of 
the piled-up PSF. These various structures to the EPIC PSF are 
discussed throughout this paper 

1.1. PSF descriptions witliin tine XMM-Newton SAS 

Historically there have been a number of descriptions of the 
EPIC PSF that have been used as part of the XMM-Newton 
Scientific Analysis System (SASU These include a number 
of two-dimensional (2-D) and one-dimensional ( 1 -D) descrip- 

' http://xmm.esa.int/sas/ 



tions. The only previous 2-D description is the 'Medium' mode 
PSF - a set of simulated images in a matrix of 6 different off- 
axis angles and 1 1 different energies. This set of 66 images is 
identical for each EPIC instrument, and there is no azimuthal 
variation included. This PSF description is the one that has 
been used for source-searching throughout the XMM-Newton 
mission, and is referred to as the 'default' PSF in this paper. 

1-D descriptions include (i) an 'Extended' mode, incor- 
porating a King profile (with core radius and slope), as a 
function of EPIC instrument, energy and off-axis angle, (ii) 
a 'High' mode, incorporating a 3-Gaussian parametrization 
of the 'Medium' mode images, and (iii) a 'Low' mode, an 
early (and now unused) single-Gaussian approximation of the 
PSF The 1-D PSFs are often used within the SAS for spec- 
tral work. These PSFs can all be found within the particular 
XMM-Newton PSF current calibration files (CCFs). These are 
named as XRT/_XPSF_«n«n.CCF, where / refers to the instru- 
ment (l=MOSl, 2=MOS2, 3=pn), and nnnn is the issue num- 
ber of the file. 

That the default (Medium) PSF is limited, in that it is the 
same for each EPIC instrument, and that no variation with 
source azimuth is included, is a drawback. This is the PSF 
that has been used for source-searching throughout the XMM- 
Newton mission, and it is likely that a number of these lim- 
itations may be responsible for many of the problems seen 
in the source-searching results and in constructing the XMM- 
Newton catalogues (e.g. 2XMM; Watson et al. 2008), the major 
problem being that a number of spurious sources are detected 
by the standard source-detection pipeline, especially close to 
bright point sources (note that post-pipeline visual screening 
does take place in the production of the XMM-Newton cata- 
logues, and sources that lie in regions where spurious sources 
are considered likely to occur, are flagged). There is therefore 
the need for a full 2-D energy- and off-axis angle-dependent 
PSF model that incorporates all aspects of the PSF structures, 
that is instrument-specific, and which accounts correctly for a 
source's azimuthal position. This paper describes this new 2-D 
PSF, and the structure of the paper is as follows: Sect. 2 de- 
scribes the construction of the PSF - the data analysis and the 
modelling. Sect. 3 describes how the resultant new PSF model 
appears in the SAS. In Sect. 4 we discuss our testing of the PSF, 
and the major findings, and in Sect. 5 we present our conclu- 
sions. 

2. Constructing the 2-D PSF: Data Analysis and 
IVIodelling 

A standard method to obtain a fully 2-D characterization of the 
PSF as a function of energy and off-axis angle across the en- 
tire field-of-view (FOV) of each of the three EPIC instruments 
would be to obtain a large number of images of appropriate 
point sources from the XMM-Newton archive, and stack these 
together in instrument/energy/off-axis angle groupings so that 
their shapes could be analysed and modelled. 

It quickly became apparent however that such an approach 
was problematic. For a specific off-axis angle, the general 
tangentially-stretched profile of the PSF that is naturally due 
to the implementation of the Wolter I optic rotates naturally 
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Fig. 1. PSFs of the 3 epic detectors (MOSl left, MOS2 centre, pn right) for the same non-piled-up source from the same obser- 
vation (2XMM J 130022. 1+282402, ObsID 0204040101, revolution 823). The images are 0.2-lOkeV, are of binning l.l"xl.l", 
and are very lightly smoothed, to accentuate the features. 




Fig. 2. (Left) A front-end view of one of the EPIC mirror modules containing the 58 co-axial mirrors shells, and the spider 
support structure used to hold the shells. (Right) The MOS2 PSF of the severely piled-up source GX 339-4 (ObsID 0204730301, 
revolution 783), showing the various PSF and other features (see text) - bad columns on the CCD are also visible as dark vertical 
lines. 



around the detector with source angle on the detector The 
PSF details that are to do with the support structure and the 
mirror deformations however - i.e. the spokes and the trian- 
gles/pentagons) do not rotate with source angle. These - re- 
ferred to collectively as the support structure features - remain 
fixed in detector angle as the source rotates around the detec- 
tor. An immediate upshot of this is that the situation is very 
complex, with every single position on each of the three EPIC 
detectors having a single unique PSF. Even for such a high- 



throughput and long-term mission as XMM-Newton, nowhere 
near enough good quality data exists to perform this full stack- 
ing analysis at every detector position. 

Consequently the new 2-D PSF system, both the modelling 
described here and the incorporation of the PSF reconstruc- 
tion into the SAS, has had to be re-designed in the follow- 
ing way: The stackings in the instrument/energy/off-axis angle 
groupings were used to construct general elliptical 'envelopes' 
(due to the optic implementation) for that particular instrument. 
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energy and off-axis angle, the support structure features (the 
spokes, triangles etc.) having been averaged and smeared out 
via the stacking procedure. 2-D spatial models were then fit to 
these stacked envelopes to obtain general envelope spatial pa- 
rameters. Then, for a particular source at a particular known az- 
imuthal angle (whether on the detector, or on the sky), the par- 
ticular azimuthal support structures required (the spokes, trian- 
gles etc.) were folded into the appropriate elliptical envelope to 
create the final 2-D PSF appropriate for that particular instru- 
ment, energy, off-axis angle, and importantly, azimuthal angle. 
This is described in more detail in the following sections. 

2.1. The Elliptical Envelope 

A major driver for this new 2-D description of the PSF was that 
it should cover the entire FOV, i.e. both on-axis and fully off- 
axis. It was necessary therefore to use data from those observ- 
ing modes that cover the full FOV - i.e. the 'full-frame' modes. 
Though smaller window modes can be useful in deriving some 
PSF parameters (see later), these modes only exist on-axis, and 
it was decided for the off-axis considerations, and for the sake 
of consistency and uniformity, to select all the data solely from 
the full-frame modes. 

Sources were selected from the 2XMM catalogue (Watson 
et al. 2008) on the basis of them: 

- Being identified within 2XMM as point-like (i.e. no exten- 
sion). 

- Having large numbers of (0.2-12keV) counts (>5000 for 
MOSl or MOS2 and >15000 for pn). 

- Having a countrate below the appropriate pile-up limit^ 
(0.70 ct/s [MOS full-frame], 6 ct/s [pn full-frame], 2 ct/s 
[pn extended full-frame]). 

- Covering the fuU range in off-axis angle. (Note that the off- 
axis angle is just the 'radial' angle relative to the optical 
axis, and any variation caused by mirror deformations (as 
opposed to support structure) with detector azimuth angle 
will be considered with the support structure features.) 

The appropriate raw data - the Observation Data Files 
(ODFs) - for each source were identified and obtained, and the 
standard SAS (v7.0) procedures ('epchain' for pn, 'emchain' 
for MOS) were run on these to create the standard calibrated 
event lists. These were then filtered for periods of high back- 
ground (solar proton flares) via standard Good Time Interval 
(GTI) files created via analysis of high-energy off-source 
lightcurves. Files (and therefore candidate sources) where a 
large amount or long durations of high background flaring was 
observed were rejected from further analysis. Large-scale im- 
ages around each source were created, and further rejections 
were made in cases of crowded fields or where chip gaps or bad 
CCD columns were seen to lie too close to the target source. 

For the sources that passed through to this stage (>250 ob- 
servations, mostly comprising useful exposures for all three 
EPICs, and several containing multiple valid sources), high- 
resolution images were constructed around each source posi- 



tion. These were constructed to be similar to the 'Medium' 
PSF images (i.e. of 1.1" resolution), and aligned such that 
the 2XMM source position lay at the centre of the image. 
These images were created for (where appropriate) each of 
the three EPIC cameras and in several different energy bands; 
0.1-l.OkeV, 1.0-2.0keV, 2.0-3.5keV, 3.5-5.0keV, 5.0-7.0keV, 
7.0-9.0keV, 9.0-11.5 keV, 11.5-15.0keV. All other sources (all 
of which were very faint and well separated from the target 
source) were excised from each image (using a circle of radius 
35"), the holes filled with the average value in the surrounding 
(40" -70") annulus. Each image was then rotated with respect 
to the source's azimuthal angle on the sky (this angle is equiv- 
alent to a combination of the position angle of the observation 
and the source's azimuthal angle on the particular detector). A 
Gaussian was fit to the rotated image to account for any tiny 
shifting in the procedure (always <<1"), and a final rebinning 
of the rotated image about this fit centre was performed. 

The final good images were collected together for each 
combination of EPIC instrument (3), energy band (8, lin- 
ear fiducial midpoints: 0.55 keV, 1.5 keV, 2.75 keV, 4.25 keV, 
6.0keV, 8.0keV, 10.25 keV, 13.25 keV) and off'axis-angle band 
(6; 0-1.5', 1.5-4.5', 4.5-7.5', 7.5-10.5', 10.5-13.5', & >13.5', 
these assigned to the fiducial values 0', 3', 6', 9', 12', 15'). 
These sets of images were then stacked together, bringing them 
together to a common reference frame. 

These stacked images were fit with a 2-D King profile using 
the beta2d model in CIAO-sherp^I : 



B{r) 



[1 + (r/ro)2]« 



where: 



r(x, y, 6) 



^[(xcose-l-ysinfl)^] + [(ycosfl - ;csin6')2]/(l - e)2 

and ro (core radius), a (index), e (ellipticity) and 9 (angle 
of ellipticity) are the model parameters. 

To the parametrisation above, a further 2-D Gaussian core 
has been added to model excess emission in the low- and 
medium-energy (E < 6 keV) PSF at the very centre of the 
MOS cameras (no such core is observed to be necessary for 
pn). The model includes the same ellipticity term through the 
definition of the r variable as above. The 2-D Gaussian function 
used {gausid model in CIAO-Sherpa) is : 



Gir) = Ae 



-41n(2)(r/FWHM)^ 



^ XMM-Newton Users Handbook : 

|http://xmni.esac.esa.int/extemal/xmm_user_suppo rt/doc umentation/uhb/| 



where FWHM is the full width at half maximum. All these 
parameters are contained within the PSF CCFs (issue num- 
bers 0011 and after). Finally, a flat background level, fixed at 
the mean value beyond a 5' radius was added to the model 
for each stacked image. For all instruments and off-axis an- 
gles, the two highest energy stacked images (at 10.25 keV and 
13.25keV) were seen to contain only very sparse data, and it 
proved impossible to obtain stable 2-D fits with sensible error 
bounds in many of these cases. Consequently, the images in 
the two highest-energy bands have been combined together for 



http ://cxc . harvard, edu/sherpa/ 
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Fig. 3. Results from a typical run of the 2-D fitting (specifically 
for pn, 6' off-axis angle, 1.5 keV). The top and middle panels 
show the data and final model respectively, at the same scal- 
ing. The bottom panel and colourbar (running from -0.2 to 0.2) 
show the residuals (data-model, the central peak of the data 
and model panels being ^5 in the colourbar units). 

all instruments and off-axis angles, and the resultant spatial fit 
parameters obtained have been applied to both energy bands 
in the resultant calibration files (at 10.25 keV and 13.25 keV). 
Future calibration, using more data, may be able to establish 
fit parameters for both energy bands separately. An example 
showing the results from a typical run of the 2-D fitting (specif- 
ically for pn, 6' off-axis angle, 1 .5 keV) is shown in Fi^ 

The actual parameter values, which are contained within 
the PSF CCFs (described in more detail in Sec l3.ll ). have 
evolved, and will continue to evolve, with iterations in the anal- 
ysis and modelling. It is possible however to discuss some 
general trends in the current parameters (CCF issue number 
0013). Dealing first with ro (core radius) and a (index), note 



that there is some degeneracy between these two parameters; 
increasing (or decreasing) both parameters simultaneously can 
sometimes lead to only very small changes in the overall pro- 
file. The two MOS PSFs behave similarly; close to on-axis, ro 
is seen to gently decrease with increasing energy, while a is 
seen to remain roughly constant with energy. The profile there- 
fore is seen to get narrower at higher energies, as expected, as 
only a fraction of the mirror shells - the smaller, inner shells - 
contribute to the higher-energy PSF. At larger off-axis angles, 
the MOS ro and a values are seen to decrease with increas- 
ing energy, but more steeply, starting from larger low-energy 
values than is seen on-axis. The pn PSF behaves slightly dif- 
ferently in that ro decreases gently with increasing energy for 
all off-axis angles, but the average ro value increases with off- 
axis angle (as expected, the PSF becoming wider off-axis). The 
pn a is seen to remain roughly constant for almost all energies 
and off-axis angles, hence, in conjunction with the gently de- 
creasing ro values with energy, the profile again gets narrower 
at higher energies. The ellipticity behaves as one would expect 
and very similarly for each instrument, getting larger with in- 
creasing off-axis angle, from on-axis to ^0.6 far off-axis. 
Also the ellipticity is seen (more prominently off-axis) to gen- 
tly increase with increasing energy in each instrument. Finally, 
for the Gaussian component (which exists only for the MOS 
cameras at £ < 6 keV), both the FWHM and the relative nor- 
malization for both MOS PSFs are seen to fall from to 6 keV, 
for all off-axis angles (i.e. to match there being no Gaussian at 
>6keV). 



2.2. The Support Structure features - The Spokes 

The radial spokes have been modelled with a flat-topped trian- 
gular function (Fig. HJ, chosen after consultations with mirror 
experts (B. Aschenbach, private communication; R. Willingale, 
private communication). 16 equally-spaced spokes have been 
included in the model. The distance between two consecutive 
peaks of the spoke function is therefore 22.5°. The shape of 
this function has been tuned in such a way that it does not 
change the integral of the radial profile (i.e., the red and the 
green areas in Fig. |4] exactly cancel out). This filtering func- 
tion is fixed in space with respect to the detector, such that any 
axis of the DETX/Y (or RAWX/Y) coordinate system lies ex- 
actly between two spokes, and is only applied to the elliptical 
envelope after the envelope has been rotated according to the 
azimuthal position of the source on the detector. 

The current spoke dimensions have been obtained from 
measuring the azimuthal angles of the spoke troughs, peaks and 
flat-tops of a number of bright sources, including GX 339-4, 
using the full energy range (0.2-12keV) and over a large ra- 
dial range (20"-100"). The current dimensions are such that 
both the half-width of the flat-top (u), and the width from base 
to peak (v) equals 0.165 times half the width of the spoke-to- 
spoke distance (as indicated in Fig|4l). It can be shown that the 
height h of the spoke above the unspoked dashed line of Fig|4] 
(if the spoke drops below the line by a depth of 1.0) is related 
to the half-width of the flat-top u and the base to peak width v 
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Fig. 4. Schematic of the flat-topped triangular function used to parameterize the PSF spokes. The blue spoking function is con- 
structed such that the red and the green areas exactly cancel out. The dashed line shows the unspoked function. The lengths 1 .00, 
0.165 (the half- width of the flat-top u) and 0.165 (the base to peak width v) are in units of half the inter-spoke distance, i.e. 
11. 25°. The figure is not to scale. 



via : 

2 

h=- 1 

2u + V 

The strength (intensity) of the primary spokes is set such 
that the spoked PSF is made up of (at most) 42% of an im- 
age formed from the flat-topped triangular function depicted 
in FigEl and the remainder is formed from an unspoked image. 
This is again based on the spatial analysis of a number of bright 
sources. The radial dependence of the intensity of the spoking 
is however not constant (see below), and this intensity value 
(42%) is the maximum value. 

16 secondary spokes, observed in very bright sources (see 
Fig|2]i, and situated half-way between the 16 main primary 
spokes, are modelled with the same spoke function, rotated by 
1 1 .25° with respect to the main spoke system, and with a (max- 
imum) intensity factor of 3.5%. 

The relative strength of the spokes is not constant with ra- 
dius but appears to vary. Owen & Ballet (2011) were able to 
characterize the radial dependence of the spokes of the MOSl, 
MOS2 and pn PSFs to large radii (up to 2') using bright, heav- 
ily piled-up sources, and using procedures developed to cor- 
rect for the eff'ects of pile-up. They found that the radial de- 
pendence of the spokes of the two MOS cameras and the pn 
camera appeared to be approximately the same. This is not sur- 
prising as the spokes and their radial dependencies are all due 
to the construction of the mirror modules, all of which were 
designed to be identical. Furthermore the spiders were con- 
structed out of Inconel, a material chosen for its thermal ex- 
pansion, which is close to that of the electrolytic nickel of the 
mirrors (de Chambure et al. 1999), thus keeping distortions to a 
minimum. This radial dependence has been modelled to match 
the Owen & Ballet (2011) on-spoke to off-spoke ratio results. 



and will be introduced into the PSF CCFs (for issue numbers 
0014 and after) as follows: For each EPIC PSF, from the cen- 
tre out to 10", there is no spoking. From 10" to 1 10" the spoke 
strength increases linearly from zero to the maximum (42%). 
From 110" to 180" the spoke strength decreases linearly from 
its maximum to zero, and there is no spoking beyond this ra- 
dius. 



Looking beyond the present model, when we compare in 
detail the azimuthal profile of good quality point source data 
with the current (CCF 0013) model (see Fig.|5]l, there are a 
number of points of interest: Though there is significant sta- 
tistical scatter, there may also be some true spoke-to-spoke 
variation. Further, it appears that there may be some varia- 
tion of the spoke model width parameters u and v with ra- 
dius - the spokes appearing to be narrower at larger radius. 
Interestingly it appears that there is a significant difiference in 
the radial dependence of the secondary spoke strength, com- 
pared to that of the primary spokes - the secondary spokes 
falling off quicker (there also appears to be significant sec- 
ondary spoke-to-spoke variation). None of these effects (nor 
their energy-dependencies) are yet contained within the cur- 
rent PSF model, but they impinge upon e.g. the on-spoke to 
off-spoke Owen & Ballet (201 1) results, and on their interpre- 
tation and modelling - note that the on-spoke to off-spoke ra- 
tio is lower at smaller radius (44"-88") than at larger radius 
(88"- 132"), partly because the secondary spokes are relatively 
stronger at smaller radius than at larger radius. It is hoped that 
these issues can be revisited in future improvements in the 
modelling of the EPIC PSFs. 



A. M. Read et al.: A New Comprehensive 2-D Model of the Point Spread Functions of the XMM-Newton EPIC Telescopes 



7 



r=44"-88" 




oL 



180 200 220 240 

Azimuth Angle (degrees) 

r=88"-132" 




180 200 220 

Azimuth Angle (degrees) 



Fig. 5. A section of the azimuthal profile (black line) from 
GX 339-4 (ObsID 0204730301) (MOSl, 0.2-12keV) in annu- 
lar extractions of (top) 44"-88"and (bottom) 88"-132". The 
current (CCF 0013) flat-topped triangular function for the pri- 
mary spokes at maximum strength (i.e. at a radius of 110") is 
shown in blue (see text). 

2.3. The Support Structure features - The Large-scale 
azimuthal modulation 

The apparent triangular shape of the MOS2 PSF can be mod- 
elled by a spatial low frequency overall modulation of the PSF 
shape. A pentagonal modulation is also present at a lower level 
in the MOS 1 and also at a very low level in the pn PSF. These 
shapes are thought to be due to distortions and irregularities in 
certain sets of mirror shells within each of the telescopes - e.g. 
the outer shells not being perfectly circular (de Chambure et al. 
1999). It is not surprising therefore that these deformities are 
different in the three EPIC PSFs. 

This modulation M(ip) is modelled as a function of the az- 
imuthal angle <p, with a multi-peaked cosine function in the 
PSFs of the cameras: 

M{4>) = A cos(N4> - 4>q) 

where A is the amplitude of the modulation, is the number of 
peaks in the modulation in a full 360° revolution, and 0o is the 
azimuthal offset of the cosine function. The current parameters 
of this function are shown in Table [T] 

The radial dependence of these large-scale azimuthal mod- 
ulations is also seen not to be constant with radius, and the 
maximum amplitude A is listed in Tabl^U The dependence is 
only roughly similar to the radial dependence of the spokes, and 
there are variations between the EPIC detectors - the pn pen- 
tagon for example only appears to be visible over quite a nar- 
row range in radius. This radial dependence of the large-scale 
azimuthal modulations (and the entire pn large-scale azimuthal 



Table 1. Current parameters of the large-scale azimuthal multi- 
peak modulation. The modulation M(<p) is modelled as a cosine 
function of the azimuthal angle 0, with A as the amplitude of 
the modulation, being the number of peaks in the modulation 
in a full 360° revolution, and (/>o the azimuthal offset of the 
cosine function, (p and <po run clockwise from north for a source 
on the sky, for an observation Position Angle, PA = (see text). 



Cameras 


A 


N 


00 


MOSl 


13% 


5 


62° 


MOS2 


45% 


3 


50° 



modulation) has not been fully calibrated yet, but will appear 
in a future issue of the PSF CCFs. 

As a final step in the construction of the 2-D PSF, a radially- 
dependent very light smoothing is applied. A smoothed image 
of the PSF is constructed using a flat boxcar of halfwidth 1.65". 
The final PSF is then a sum of / times this smoothed image plus 
(1 - /) times the original unsmoothed image, where / varies 
linearly from at the source center to unity at 8.8". Though this 
does not currently directly match the more elliptical nature of 
the PSF at larger off-axis angle, this will be investigated further, 
and the smoothing effect is very small, has been applied merely 
to de-pixelate the spoke structures at large radius, and does not 
have any significant large-scale effects (e.g. on radial profiles). 

2.4. Additional features 

The new 2-D PSF system has been designed to be modular, and 
it is quite straightforward to add in further complexities. There 
are a number of additional features to the EPIC PSFs that will 
be added, it is intended, to future updates of the system (some 
may just require updates to the CCFs, others may require soft- 
ware and infrastructure changes). The radial dependencies of 
the large-scale azimuthal modulations, as discussed above, is 
one. It is very probable, as the large-scale azimuthal modula- 
tions are due to deformations only in certain mirror shells, that 
there are energy-dependencies to this effect, and these will also 
need to be incorporated. A further task is to include the effects 
of the dark lanes visible e.g. in Fig|2] due to the electron deflec- 
tor, mounted after the rear aperture of the mirrors and whose 
legs align with those of the front-end spider It should be pos- 
sible, once calibrated, to insert this effect into the system, via a 
similar azimuthal filtering technique as for the spokes. Features 
due to OOTs have not been considered to be part of the PSF. 
Though instrument-dependent, they are also mode-dependent, 
and in terms of the source-detection software, they are mod- 
elled via an alternative route (at least for the pn, where the OOT 
signal is much stronger). The OOT features here have, in a sim- 
ilar manner to the azimuthal structures (the spokes, triangles 
etc.), been diluted and smeared out via the stacking procedure 
prior to the envelope fitting. Furthermore there is set within the 
PSF system the possibility to allow for any azimuthal variations 
in the envelope 2-D King and Gaussian parameters. Currently 
the CCFs are set such that these parameters do not vary with 
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source azimuth, but future calibration may require changes, 
e.g. to account for variations due to the obscuring RGAs in the 
EPIC-MOS cameras. 

3. Results: The 2-D PSF in the XMM-Newton SAS 

The new 2-D PSF is contained within all XMM-Newton 
SAS PSF Current Calibration Files (CCFs) of the form 
XRT/_XPSF_n«n«.CCF, where nnnn, the issue number of the 
file, is 0011 or higher. The relevant parameters are contained 
within the ELLBETA_P ARAMS extension. 

A number of elements are used to describe the new 2-D 
PSF in these CCFs, as summarized here: 

- a (currently) azimuthal-angle-independent parameteriza- 
tion of the elliptical PSF envelope as a function of instru- 
ment, energy, and off-axis angle. An additional Gaussian 
'core' is added to the PSF envelope of the MOS cameras 
for low and medium energies (E < 6 keV) only. 

- an azimuthal filter of the elliptical envelope, which de- 
scribes the effect of the spokes created by the spider sup- 
porting the 58 co-axial mirrors of each telescope. Primary 
and secondary spokes are included. The radial dependency 
of the strength of these spokes is to be included in version 
0013 of the CCFs. The parameters of this filtering are cur- 
rently hard-coded into the CAL software, but will appear 
as FITS keywords in version 0013 of the CCFs. The filter- 
ing to account for the dark lanes will appear in future CCF 
releases. 

- a further azimuthal dependency of the overall PSF enve- 
lope, which is responsible for the apparently triangular 
shape of the M0S2 PSF. Lower level similar effects are 
present in the PSFs of MOS 1 and pn also. Currently, this 
con-ection is applied to the PSFs of MOS 1 and MOS2 only. 
The parameters of this correction are currently hard-coded 
into the CAL package, but will appear as FITS keywords 
in version 0013 of the CCFs. The pn correction and the ra- 
dial dependencies for all three will be included in future 
releases of the CCFs. 

A scheme of the different physical ingredients of the 2-D 
PSF parameterization and how they are combined to produce 
the final PSF is shown in Fig.|6] For a comparison with an ac- 
tual EPIC off-axis source, Fig.|7]shows a very bright, ^ 4' off- 
axis angle MOS2 point source and the equivalent PSF model 
at a similar off-axis and the appropriate source azimuthal posi- 
tion. 

Once the 2-D PSF was incorporated into both the CCF sys- 
tem and the SAS, a next step was to compare it with the default 
PSF in terms of how well the SAS-constructed PSFs model the 
true source structures. To this end, full SAS calibrated event 
list creation, image formation and detection chain analysis (as 
per the standard pipeline) was performed on a large number of 
fields containing bright, non-piled-up point sources, detected 
by EPIC, and (where available and where observing mode con- 
straints allowed) in all three EPIC cameras. The detection chain 
analysis was performed once using the default (MEDIUM) PSF 
and once using the 2-D (ELLBETA) PSF. Next, taking great 



care to co-align each source image and each detection-chain- 
produced model PSF image, then for each source, both the 
source and the model counts within each of 30 radial bins (each 
4" wide) and 32 azimuthal bins (each 1 1.25° wide, arranged to 
be alternately fully on-spoke, then fully off-spoke) were com- 
puted. The co-alignment ensured that each azimuthal bin cor- 
responded to the same azimuth on the PSF. This procedure was 
performed for 80 bright, non-piled up point sources, and the 
values were stacked together, in several energy bands and off- 
axis angle groupings. All this analysis was performed using 
both the 2-D PSF and the default PSF 

FiglS] shows the results of this analysis, for the example of 
M0S2, on-axis, 0.5-1 keV, arranged in the circular grid used 
of 30 radial bins combined with 32 azimuthal bins. The top row 
corresponds to the 2-D PSF, and shows (left to right): (1) the 
stacked counts in the data, (2) the stacked counts in the model, 
(3) a ;t'^-like statistic, calculated from the stacked data counts 
and stacked model counts as : 

, +(stacked data - stacked model)^ 

X — ' 

stacked model 

and (4) r, a measure of whether the model consistently under- 
or overestimates the data, and calculated as the sum of the in- 
dividual data-model 'signs' via : 




data - model 
|data - modell ' 



which ranges between -i-l (white, always underestimates) and 
-1 (black, always overestimates). The bottom row of Fig|8] 
shows the equivalent plots for the default PSF (and hence the 
two stacked data plots are identical). The colourbar (-20 to -h20) 
corresponds to the plot. 

A number of conclusions can be drawn from this figure. 
Similar results are seen for MOSl and pn, though some (e.g. 
to do with the large-scale azimuthal filtering) are harder to 
see. Very obvious is that the default PSF, due to its limita- 
tions discussed in Sec ll.ll performs very poorly with regard 
to the spokes; the on-spoke regions are very underestimated 
by the default model, and the off-spoke regions are overesti- 
mated. Secondly, the default PSF performs extremely badly in 
modelling the large-scale azimuthal features - here, the MOS2 
triangle (note the very large triple-peaked discrepancies in the 
and r plots for the default PSF). The 2-D PSF however, per- 
forms very much better Much less variation inx^ and r is seen 
on- and off-spoke. The plot does suggest though that the 2-D 
PSF relative spoke strength may be too strong. The most up-to- 
date SAS and the latest PSF CCF files that are available for this 
analysis (version 0012) however, do not yet include the radial 
dependency of the spoke strength, and the future inclusion of 
this will act to decrease this effect, and improve further the 2-D 
PSF modelling. The 2-D PSF modelling of the MOS2 triangle 
is also very much improved over the default PSF. Further re- 
finements to this component include the modelling of both its 
radial-dependence and its energy-dependence. 

3.1. Description of tine PSF CCF parameters 

The FITS extension ELLBETA_PARAMS in the PSF CCF 
XRT/_XPSF_0011.CCF (and later) files contains the parame- 
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Fig. 6. The eight main steps in the formation of the full 2-D PSF for a source in a given instrument, of a given energy and at a 
given off-axis and azimuthal angle: The King (beta2d) component [1] is constructed, then the Gaussian (gaus2d) core [2] is 
constructed, and these are added [3] in the correct ratio (the CCF parameters in steps 1-3 are all functions of instrument, energy 
and off-axis angle). Then this is rotated [4] according to the azimuthal position of the source on the detector, and only then 
are the radially-dependent primary [5] and secondary [6] spoke structures azimuthally filtered in, using a flat-topped triangular 
function. Finally, the large-scale azimuthal modulation (a function of EPIC instrument) is filtered in [7], and the very light 
radially-dependent smoothing applied [8]. The example shown is for MOS2, at an energy of 1.5 keV, an off-axis angle of 9', and 
a source azimuthal position of 30°. 



ters describing the elliptical envelope and the Gaussian core. 
This extension contains four columns: 

- ENERGY: the energy (in eV) to which the parameters refer 

- THETA: the off-axis angle (in radians) to which the param- 
eters refer 

- PHI: the azimuthal angle (in radians) to which the parame- 
ters refer 

- PARAMS: an array, containing the parameters of the 2-D 
King (beta2d) plus Gaussian (gaus2d) function: 

- 1: the King core radius (ro), in arcseconds 

- 2: the King power-law slope (a) 

- 3: the ellipticity (e) (of both the King and the Gaussian 
components) 

- 4: the Gaussian Full Width Half Maximum (FWHM), 
in arcseconds 

- 5: the normalisation ratio of the Gaussian peak to the 
King peak (A^) 



4. Discussion: Testing of the PSF 

Though this paper is primarily concerned with a description of 
the new 2-D PSF and its construction etc., it is very useful here 
to discuss some of the testing, especially as regards the major 
reason why the new 2-D PSF was originally desirable - i.e. to 




Fig. 7. A very bright, slightly piled-up, x 4' off-axis angle 
MOS2 point source and the equivalent PSF model at a simi- 
lar off-axis and the appropriate source azimuthal position. 



see if it could reduce the number of spurious sources detected 
by the standard source-detection pipeline, where the relatively 
poor description of the PSF can lead to large residuals in the 
constructed data-minus-model images, leading in turn to the 
detection of spurious sources. Also of major interest are the 
astrometry issues that have arisen as a consequence of these 
first tests. These issues are discussed in the following sections. 




Fig. 8. Comparing the SAS-produced 2-D PSFs and default PSFs with true source data. Left to right: (1) the stacked counts in the 
data, (2) the stacked counts in the model, (3) a;t-^-like statistic, (4) r, a measure of whether the model consistently (white) under- 
or (black) overestimates the data (see text for all details), for (top) the 2-D PSF and (bottom) the default PSF. The colourbar 
corresponds to the;^^^ plot. The 30 radial bins are each 4" wide, and the 32 azimuthal bins are each 11. 25° wide, and arranged to 
be alternately fully on-spoke, then fully off-spoke. The example shown is for MOS2, on-axis, 0.5-1 keV. 



4.1. Source-searching and spurious sources 

A comparison of the performance of the new 2-D PSF 
('Ellbeta') with that of the current existing default PSF 
('Medium') has been performed, the primary objective be- 
ing to establish whether the 2-D PSF better suppresses the 
detection of spurious sources. The analysis first focussed on 
a set of 108 'problem' observations containing examples of 
various types of situations where spurious sources had pre- 
viously been found to be present. These were mainly chosen 
to be around non-piled-up bright sources, but there were also 
several fields chosen involving bright piled-up sources, off- 
axis sources, extended sources, out-of-time (OOT) events from 
pUed-up sources and single reflection arcs. A separate 'ran- 
dom' set of 83 observations was also processed as a reference, 
as was a 'clean' subset of 40 fields, free of any 'problems', 
extracted from the 'random' subset. 

Standard calibrated event file creation, and background 
cleaning was performed on the datasets, and standard images 
were extracted from the data in the 5 usual energy bands used 
for EPIC source-detection in the 2XMM catalogue (Watson et 
al. 2008). The main source-detection task emldetect (v5.15) 
was then used, once using the default PSF, and once using the 
2-D PSF. Importantly, exactly the same input files (specifically 
the source images, background images, exposure maps, detec- 
tor masks and eboxdetect input lists) were used in the default 
and 2-D runs, ensuring that no changes were introduced due to 
e.g. spatial randomization of the events, and that any changes 
seen would be due solely to the PSF usage. 



Table 2. Total numbers of emldetect detections with the 2-D 
and default PSFs, and the percentage change in numbers, for 
various 'problem' cases, and for the 'random' and the 'clean' 
samples (see text). Negative percentage changes indicate that 
the 2-D PSF performs better at avoiding likely spurious sources 
in problem fields. 



Set 


Problem 


N 


N 


(2-D - def)/def 




type 


(2-D) 


(def) 


perc. change 


Problem 


non-piled-up 


4286 


4491 


-4.6% 




piled-up 


1171 


1255 


-6.7% 




off-axis sources 


1372 


1388 


-1.2% 




OOTs (piled-up) 


132 


135 


-2.2% 




extended sources 


896 


893 


-hO.3% 




reflection arcs 


1194 


1172 


+ 1.9% 


Problem 


Total 


9417 


9700 


-2.9% 


Random 


Total 


6395 


6524 


-2.0% 


Clean 


Total 


2570 


2575 


-0.2% 



The numbers of emldetect detections obtained using the de- 
fault and the 2-D PSFs for the 'problem' fields (split into the 
various problem groups), and for the 'random' fields and the 
'clean' fields are shown in Table|2] 

Overall for the problem cases, the number of detections 
with the 2-D PSF is lower (by 2.9%) than for the default PSF 
This is dominated by the large reductions in the on-axis non- 
piled-up and piled-up cases. The same analysis applied to the 
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random cases (where many similar 'problems' still exist) yields 
a smaller reduction than for the problem cases, while for the 
clean cases, there is almost no change. This all supports the 
idea that the 2-D PSF performs better at avoiding likely spuri- 
ous sources in problem fields. In fact, if one removes the clean 
subset from the random set, the excess of detections by the de- 
fault PSF in the remaining 'non-clean' cases is 3. 1 %, consistent 
with the 'problem' 2.9% in Table|2] Two examples showing the 
improvement in the source detection around bright sources are 
shown in Fig|9] 

We then focussed on the detection statistics around the (typ- 
ically brighter) objects that give rise to the excessive source 
detections. Using the set of problem fields, source counts were 
recorded for 3 circular apertures centred on the problem source 
of interest, with radii of 1', 2' and a radius determined by eye 
that enclosed the visible spoke structure. For insignificantly 
piled-up sources (examining 53 fields, with essentially one 
bright object per field) we find that use of the 2-D PSF reduces 
the number of detections around bright sources by a;23% on av- 
erage, and is not noticably dependent on the aperture. Applying 
the same analysis to significantly piled up sources (14 sources) 
and to bright off-axis sources (21 sources) yields changes of 
between 26% and 31%. 

To examine more closely the distribution of sources around 
problem objects, a visual analysis was conducted on a subset 
of objects, counting sources that were either associated with 
the spoke features, or suspect objects found elsewhere in the 
2' aperture but generally close to the core of the source. The 
primary source detection itself was excluded. This is inevitably 
a subjective analysis and does not take account of the depen- 
dence on detection likelihood, but it provides a useful assess- 
ment of the impact of the 2-D PSF. We find that for sources 
detected on spokes, overall the 2-D PSF finds ;J,42% fewer de- 
tections than the default PSF (or >92% fewer if we consider 
only detections on spokes that are not common to both PSFs). 
Likewise, for sources near the core of the primary source, the 
2-D PSF yields at least 22% fewer detections, or >69% fewer 
if we exclude the common detections. 

Thus, we conclude that use of the 2-D PSF substantially 
reduces the number of detections compared to the default PSF 
and that this reduction is primarily due to fewer detections on 
the spokes and in the core regions of bright sources. These are 
precisely the components of the PSF which are being better 
modelled by the 2-D PSF. These model improvements suppress 
data-minus-model residuals that have hitherto given rise to spu- 
rious detections when using the default PSF. 

4.2. Astrometry issues 

A major outstanding concern in the comparison of the 2-D 
PSF output with the default PSF output was that large and 
significant positional shifts are observed. Mean differences of 
x!+0.8"in RA and x!-0.8"in Dec (2-D values minus default val- 
ues) are observed between the 2-D and default PSF results on 
any large set of observations. Fig[TO]shows the distributions of 
the RA and Dec (2-D-default) offsets for a sample of 70 ob- 



servations (incorporating and extending on the clean sample of 
the previous section). 

It was not known whether this shift was due to the 2-D 
PSF or to the default PSF (or to both). One source of con- 
fusion and possible error is the fact that the default PSF (as 
described briefly in Sect. 1.1) consists of a set of images of di- 
mensions 512x512 and pixel size l.l"xl.l". It was noted that 
the assumption of where the centre of one of these images is, 
and the propagation, correctly or incorrectly, of this assumption 
through the relevant SAS subsystems (PSF and image genera- 
tion, PSF rotation, source-searching etc.) could lead to possible 
errors. This thinking also indicated that the rotation of the PSF 
could be vital in pinning down the root of the positional off- 
set problem, and indeed, when the (2-D - default) changes in 
RA and Dec are plotted against the source angle on the sky 
(FigfTTll. the H-0.8" and -0.8" offsets are resolved into not only 
offset, but sinusoidal variations, indicating that (i) not only is 
there an offset problem between the two PSFs, but (ii) the ro- 
tation of one (or both) of the PSFs contains an additional sys- 
tematic sinusoidal error. 

In order to evaluate which of the PSFs is at fault, it is nec- 
essary to cross-correlate the X-ray positions obtained with the 
two PSFs with good-quality optical positions. Here we need 
to introduce the boresight matrix. This was calculated from 
fields of many bright sources, e.g. OMC2/3 (Kirsch 2004), and 
subsequently refined using 430 individual fields (Altieri 2004), 
cross-correlating the source X-ray positions (note, obtained us- 
ing the default PSF) with their optical positions, and observing 
which particular offset in RA, Dec and position angle is seen 
to significantly optimize the correlation. This final calibrated 
boresight misalignment matrix has been applied to the determi- 
nation of the X-ray sky position of every X-ray event in every 
calibrated event file. 

It is therefore the case that any translation or offset problem 
that is due to the default (Medium) PSF is 'corrected for' and 
calibrated out when implementing the current boresight mis- 
alignment matrix. As such, it is expected that the X-ray-optical 
source positional offsets for the default PSF will be centred 
around zero, and indeed, this is seen to be the case. The 2-D 
PSF offsets are furthermore seen to be centred around the ob- 
served (FigslTO]&[II]l -1-0.8" and -0.8" off^sets. 

The clean samples of source detections (using both the 2-D 
PSF and the default PSF) were cross-correlated with the Sloan 
Digital Sky Survey Quasar Catalog (Schneider et al. 2007) 
to a match radius of 10". In almost every case where there 
was a match, it was the only match. We looked at the cross- 
correlations using the original (emldetect-obtained) X-ray po- 
sitions and using the X-ray positions obtained by the SAS pro- 
gram eposcorr. This program, eposcorr, performs a similar task 
to the boresight misalignment matrix calculation in that, for 
each observation, it correlates the input (emldetect) source po- 
sitions with the positions from optical source catalogues, and 
checks whether there are offsets in RA, Dec and position angle 
which optimize the correlation. If there are optimum offsets, 
these are then used to correct the input source positions which 
are then added as separate columns to the input X-ray source 
list. 
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Fig. 9. Example output from a full all-EPIC source-detection analysis of ObsIDs 0107660201 (left) and 0302850201 (right). The 
blue circles show the sources detected using the default PSF and the yellow circles show the sources detected using the 2-D PSE 
Many spurious sources previously detected by the default PSF in the spokes of the central bright source are now not detected by 
the 2-D PSE 
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Fig. 10. Distributions of the (2D-default) change in returned (left) RA and (right) Dec of detected sources using the 2D (Ellbeta) 
and default (Medium) PSE for a sample of 70 clean observations. 



When the (X - QSO) differences in RA and Dec be- 
tween the X-ray positions and the optical quasar positions are 
plotted against source angle on the sky for both the default 
PSE (EigHJli and the 2-D PSE (EigHS] both figures using the 
eposcorr-corrected output), several things are evident. 

Firstly, regarding the off'set problem, though the corre- 
sponding figures prior to eposcorr (not shown) do indeed show 
the default PSE distribution to be centred around zero offset, 
and the 2-D PSE distribution to be centred around the -1-0.8" 
and -0.8" offsets, the 2-D PSE eposcorr output (EigfT3]l is cen- 
tred around zero - i.e. when it runs successfully, eposcorr is 
able to correct the off'set problem. Secondly, regarding the si- 



nusoidal problem, it is observed that it is for the default PSE 
case (EigfTSll that the large sinusoidal variations are seen. Eor 
the 2-D PSE (EigfTSll the variations are very much smaller. 

Looking in more detail, it was suggested earlier that the si- 
nusoidal problem could be related to the default PSE CCE im- 
ages, and the assumptions regarding where the centres of these 
images are, and how these assumptions are propagated through 
the entire PSE generation and source-searching system. This 
now does appear to be the case, and added to the lower pan- 
els of EiglT2]are simple calculated cosine curves of what one 
would expect to see, given the simple off'set of half the diago- 
nal of a single default PSE image pixel, rotated around the 360 
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Fig. 11. (2D-default) offsets in (left) RA and (right) Dec of detected sources using the 2D (Ellbeta) and default (Medium) PSF 
plotted against the source angle on the sky (anti-clockwise from north) for the sample of 70 clean observations as in FigfTOl 
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-4 -2 2 4 -4 -2 2 ^ 

Default PSF (X-QSO) dRA Default PSF (X-QSO) dDec 



Fig. 12. (Top) X-ray position - QSO optical position offsets in (left) RA and (right) Dec of default (Medium) PSF-detected 
X-ray sources with a QSO match, plotted against the source angle on the sky (anti-clockwise from north). X-ray positions are 
after eposcorr (see text). (Bottom) The same data, binned into 25 equi-angular bins and averaged. 



degrees of the detector The data is seen to match this simple 
model extremely well. It is therefore undoubtedly the default 
PSF that is the cause of the sinusoid effect, and this problem has 
existed for the entirety of the XMM-Newton mission. The full 
positional capabilities of EPIC therefore have not been used 



to this date. The improvement here is such that, selecting the 
67% of sources with a very good X-ray positional error (cen- 
troid error <1"), the mean X-ray-QSO positional offset is re- 
duced from 1.13" (default) to 0.94" (2-D), and the percentage 
of these sources with an X-ray-QSO positional offset less than 
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Fig. 13. (Top) X-ray position - QSO optical position offsets in (left) RA and (right) Dec of 2-D (Ellbeta) PSF-detected X-ray 
sources with a QSO match, plotted against the source angle on the sky (anti-clockwise from north). X-ray positions are after 
eposcorr (see text). (Bottom) The same data, binned into 25 equi-angular bins and averaged. 



1" increases from 52% (default) to 68% (2-D). The percent- 
age of cases where it is the 2-D PSF X-ray position that is the 
closest to the QSO position is 70%. 

Although the large sinusoids present in FiglT2]have been 
largely removed in Fig[T3] they do not appear to have been 
completely removed (and are at the 0.2"-0.3" level). One 
should bear in mind that, if there is any inaccuracy at all present 
in the true positioning of the PSF centre, and where it is as- 
sumed to be in the PSF system, then a very similar sinusoidal 
variation will be seen with source angle on the sky. The fact 
that some residual curvature possibly exists in Fig[T3]may in- 
dicate that the situation is not quite perfect, and there may still 
be a small residual misalignment of the centering of the PSF 
system and the true PSF system. 

At present, the 2-D PSF only returns the optimum positions 
when eposcorr is run. Eposcorr only provides gross field shifts 
and cannot correct for the residual sinusoidal effects. Also, the 
eposcorr task can not be run on every EPIC dataset, for instance 
in cases where there are very few (or very faint) X-ray sources 
in the field. In order for the full positional improvements of the 
2-D PSF to be usable therefore, a revised boresight misalign- 
ment matrix first needs to be calibrated, tested and incorpo- 
rated correctly into the SAS (related SAS changes may also be 
required). This work, building on the above tests, is currently 
underway, and a future version of the SAS will include these 
improvements in the positional accuracy of EPIC. 



5. Conclusions 

A new and fully comprehensive fuU-field-of-view (FOV) 2-D 
model of the point spread functions (PSFs) of the three XMM- 
Newton EPIC telescopes has been constructed. This has been 
performed via the stacking, and bringing-together to a com- 
mon reference frame, of a large number of good quality, long- 
exposure, non piled-up, bright point sources from different po- 
sitions within the full FOV of each EPIC detector The resultant 
general PSF envelopes were then azimuthally filtered to con- 
struct the primary and secondary spoke structures (plus their 
radial dependencies) and the large-scale gross azimuthal PSF 
deformations that are observed. The PSF model also includes 
an additional Gaussian core, which accounts for (at most) 2% 
of the enclosed energy flux in the EPIC-MOS cameras. 

This PSF model is available for use within the XMM- 
Newton Science Analysis System (SAS) via the usage of 
Current Calibration Files (CCFs) XRT/_XPSF_0011.CCF and 
later versions. The modular nature of the new PSF system al- 
lows for further corrections and refinements in the future. 

Initial EPIC source-searching tests using this new PSF 
model indicate that it performs significantly better with regard 
to the major problem with the previous PSF; that of large num- 
bers of spurious source being detected in the wings of, or close 
to bright point sources. The numbers of these spurious sources 
detected with the new PSF model are greatly reduced. 

These tests also uncovered a systematic error in the previ- 
ous PSF system, such that returned source RA and Dec val- 
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ues were seen to vary sinusoidally about the true position 
with source azimuthal position angle on the sky. This error in 
the previous PSF system (the amphtude of the sinusoid being 
ssO.S" in RA and in Dec) has existed since the beginning of the 
XMM-Newton mission, and affects all the EPIC positional de- 
terminations performed thus far. Usage of the new PSF results 
in a much smaller amplitude sinusoid, and therefore an im- 
proved positional accuracy. SAS changes, including a revised 
boresight misahgnment matrix (presently under construction) 
are required to make fuU use of the improvements in the posi- 
tional accuracy of EPIC introduced by the new PSF. 
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